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INTRODUCTION  AND  BRIEF  HISTORY 


The  original  work  on  microwave  waveguide/cavity  discontinuities  in  the 
electromagnetic  field  sense  grew  out  of  applications  of  variational  principles 
(formerly  characteristic  of  quantum  mechanics)  to  electromagnetic  theory. 
This  work  began  during  World  War  II  at  the  Massachusetts  Institute  of 
Technology  Radiation  Laboratories;  Marcuvitz  (Reference  1)  and  Schwinger 
(Reference  2)  were  responsible  for  the  majority  of  theoretical  developments. 
They  exploited  the  relationship  between  electromagnetic  boundary  value 
problems  and  equivalent  microwave  networks.  Theoretical  determination  of 
lumped  circuit  parameters  was  emphasized  that  employed  variational  methods, 
integral  equation  methods,  equivalent  static  techniques  using  conformal 
mapping  (Reference  3),  and  the  transform  or  Wiener-Hopf  technique. 

All  of  the  above  techniques  assume  that  the  terminal  plane  is  far  enough 
away  from  the  discontinuity  to  consider  propagating  mode  behavior  only — the 
evanescent  modes  are  ignored.  If  we  desire  accurate  values  for  the  reflection 
and  transmission  coefficients  of  the  dominant  mode  close  to  the  Junction  or 
discontinuity  plane  (or  for  two  or  more  closely  spaced  discontinuities),  any 
analysis  must  include  the  evanescent  modes  explicitly. 

With  the  advent  of  large  digital  computers,  explicit  inclusion  of 
evanescent  modes  became  possible,  and  the  first  attempts  using  mode¬ 
matching  were  reported  (References  4  through  7).  It  was  assumed  that  for 
any  junction/discontinuily  in  a  waveguide  where  the  geometrical  shape 
resulted  in  separable  solutions  of  the  Helmholtz  equation,  the  fields  in  each 
region  could  be  represented  by  infinite  eigcnscries,  which  were  orthogonal 
and  normalized.  Using  the  continuity  of  tangential  electric  and  magnetic 
fields  over  apertures  along  with  n  x  E  =  0,  n  •  B  =  0  for  perfectly  conducting 
obstacles,  the  boundary  conditions  are  formed  in  terms  of  products  of 
eigenfunctions  and  associated  constant  coefficients,  which  are  the  modal 
reflection  and  transmission  coefficients.  These  coefficients  arc  evaluated  by 
invoking  orthogonality  properties  of  the  eigenfunctions  and  integrating  over 
appropriate  parts  of  the  discontinuity  plane. 

This  approach  was  used  by  many  investigators,  not  only  for  metallic  but 
also  for  surface-wave  waveguide  discontinuities  (Reference  8).  However,  two 
major  problems  exist: 

1.  The  original  scries  that  satisfies  both  the  differential  equation  and  the 
boundary  conditions  is  infinite,  and  it  must  remain  so  if  an  exact  answer  is 
required.  Under  a  very  few  special  circumstances  (Reference  7),  it  is  possible 
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to  invert  the  resulting  infinite  dimensional  matrix  to  solve  the  problem; 

however,  in  general,  this  is  impossible.  Thus,  to  obtain  numerical  values  of 
the  coefficients,  the  series  must  be  truncated  at  some  point.  No  one  seemed  to 
know  the  answer  to  the  question:  "How  many  modes  give  a  sufficiently 
accurate  result?"  Finally,  Lee,  Jones,  and  Campbell  (Reference  9)  shed  light 
on  this  in  a  .  very  careful  treatment  of  the  eigenfunction  form  of  mode- 
matching  and  its  relationship  to  an  integral  equation  method.  They  found  that 
for  certain  types  of  waveguide  bifurcation  problems,  the  lack  of  convergence 
in  the  solutions  was  caused  by  the  particular  form  of  mode-matching  used  and 
could  be  removed  by  using  the  moment  method.  However,  in  the  class  of  iris 
discontinuity  problems,  they  found  difficulties  with  convergence  even  when 
using  the  moment  method. 

2.  Once  a  truncation  value  for  the  series  has  been  chosen,  the  ratio 
between  the  number  of  modes  in  different  regions  must  be  considered.  This 
led  to  the  phenomenon  of  relative  convergence  (References  4  and  7)  in  which 
the  coefficients  converge  to  different  values  depending  upon  the  number  of 
modes  taken  in  each  region  and  the  ratio  of  the  modes  between  regions. 

The  failure  to  satisfy  the  edge  condition  along  the  discontinuity  plane  was 
thought  to  cause  relative  convergence  (References  10  and  11).  It  has  been 

claimed  by  many  that  there  is  a  lack  of  uniqueness  in  the  solutions  found  from 
using  only  continuous  tangential  fields  across  a  junction.  Only  when  the  edge 
condition  is  added  are  we  guaranteed  to  find  the  one  modal  ratio  that  gives  the 

physically  correct  result.  However,  as  far  as  we  know,  this  never  has  been 

proved  mathematically;  several  investigators  have  claimed  that  the  modal  ratio 
that  causes  the  system  of  linear  equations  to  be  most  "well-conditioned"  leads 
to  the  correct  result  (Reference  12). 

The  least-squares  boundary  residual  method  (LSBRM)  was  proposed  to 
avoid  completely  the  problem  of  relative  convergence  (References  13  through 
15).  Via  numerical  experimentation,  we  will  show  that  this  method,  when 
combined  with.^  traditional  mode-matching,  converges  to  the  same  answer 
regardless  of  the  modal  ratio  assumed.  This  method  not  only  avoids  the 
convergence  problem  but  has  several  other  advantages. 

1.  The  reflection  and  transmission  coefficients  can  be  found  separately, 
thus  cutting  the  size  of  the  matrix  to  be  inverted  by  half  (through 
partitioning). 

2.  The  matrices  to  be  inverted  arc  Hermitian,  thus  providing  further 
numerical  savings. 

3.  The  technique  is  characterized  by  an  error  function,  which  must  be  a 
minimum  when  the  correct  result  is  obtained. 

The  LSBRM  has  one  drawback — it  converges  slowly.  Attempts  to  include  a 
"convergence  factor"  to  speed  convergence  have  not  been  successful  in  the 
sense  that  no  a  priori  way  of  calculating  this  factor  is  known.  However,  slow 


4 


NWC  TP  6941 


convergence  is  a  small  price  to  pay  for  a  method  in  which  a  very  stable 
scattering  matrix  can  be  formulated  that  satisfies  the  traditional  boundary 
conditions  in  a  least-squares  sense. 


We  have  used  this  technique  on  two  canonical  problems:  (1)  the  infinitely 
thin  capacitive  iris  in  a  parallel-plate  waveguide  and  (2)  the  H-plane  parallel- 
plate  waveguide  step. 

Excellent  results  have  been  obtained  for  both.  In  the  following  section, 
the  general  theory  of  the  LSBRM  will  be  discussed  and  a  scattering  matrix 
formulation  (in  the  least-squares  sense)  will  be  developed.  The  above 
problems  will  be  discussed  in  the  Canonical  Problems  section,  and  results  will 
be  shown.  (The  mathematical  details  of  these  two  problems  are  in  Appendixes 
A  and  B,  respectively.)  Prior  to  our  conclusions,  we  will  indicate  how  the 
relative  convergence  phenomenon  is  avoided  by  the  LSBRM. 


GENERAL  THEORY 


In  general,  we  will  consider  a  junction  between  two  waveguides  (cavities) 
in  the  following  way  (see  Figure  1).  The  a(‘^’s  (i  =  1,2)  are  known  inputs,  and 
the  b(‘)’s  (i  =  1,2)  are  unknown  outputs  for  the  junction  at  z  =  0.  Using  the 
scattering  matrix  representation. 


b  =  S  a 


(1) 


where  b  and  a  arc  complex  column  vectors,  and  S  is  a  complex  matrix,  and 


'Sn  :Si2' 
1 

f 

<  S21  ]  S22  y 

in  partitioned  form. 

The  object  of  the  analysis  is  as  follows.  We  wish  to  form  a  scattering 
matrix  [based  on  mode-matching  and  least-squares  satisfaction  of  the 
boundary  conditions  at  the  discontinuity  (z  =  0)]  that  docs  not  exhibit  the 
relative  convergence  phenomenon.  This  is  done  to  solve  for  the  unknown 
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outputs.  We  begin  with  the  usual  infinite  sums  of  eigenfunctions  as  tangential 
electric  fields,  i.e., 


£(1) 

Man 


(1)  e-y™  "  + 

m  ''  ^ 


(1) 


] 


.(1) 

■'m 


<  0 


m=l 


(3a) 


.(2)_ 


+  a 


(2) 


■Py. 

e^p 


.(2) 


z  >  0 


p=i 


(3b) 


and  magnetic  fields 


m=l 


(4a) 


H 


(2) 

tan 


.(2). 


-  a 


(2) 


h 


(2) 


>  0 


p=i 


(4b) 


where  c(‘^  and  h(‘)  (i  =  1,2)  are  the  electric  and  magnetic  eigenfunctions, 
respectively,  and  y(‘)  (i  =  1,2)  is  the  modal  propagation  constant  along  the  z- 
axis  (in  each  region  of  the  guiding  structure);  m  and  p  represent  the  number 
of  modes  to  be  used  in  regions  1  and  2,  respectively. 

The  tangential  boundary  conditions  at  z  =  0  are 


Elan  ~  Etan  ;0<X<b  (5a) 

;  0<x<b  (5b) 

E,'”  =  0  ;  b<x<a  (5^, 
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Wc  must  recall  also  that 


hJ”  =  J,  on  z  =  0 


,  b  <  X  <  a 


but  it  is  not  necessary  to  include  this  condition  explicitly  in  order  to  obtain  a 
unique  solution  to  Equation  1.  This  and  the  associated  normal  boundary 
conditions  will  be  discussed  in  the  Canonical  Problems  section  and  Appendix  C. 

Using  Equations  3  and  4  in  5,  the  boundary  conditions  become 

09  «0 

+  Dm  J  Cm  =  ^  [bp  +  hp  J  Cp  ;  0  <  X  <  b 

m=l  p=l  (6a) 


Xta!^>-b<”lh”>.Xtb“>-aflh^^>  :  0<x<b 

m=l  P=1  (6b) 


S  =  0  ;  b  <  X  <  a 

*"=•  (6c) 


Arranging  the  b's  on  one  side  and  the  a's  on  the  other,  Equations  6  may  be 
written  in  matrix  form  as 


,(1) 


'm 


-h 


(1) 

m 


-e 


-h 


(2) 


(2) 


V 


0 


Dm 

y 

h(2) 

Wp  J 

-e 


-h 


(1) 

m 

(1) 

m 


.(1) 

'm 


e 

-h 

0 


(2) 


P 

(2) 

P 


y 


y 


(7) 


This  leads  to  a  matrix  equation  of  the  abstract  form 


L  b  =  M  a 
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Many  investigators  (References  5  through  7)  simply  assumed  that  by  applying 
the  orthogonality  properties  of  the  eigenfunctions  and  integrating  all 
boundary  conditions  over  the  appropriate  cross  section  at  z  =  0  the  unknown 
outputs  are  given  by 


L'^  Mdc)  a 


(9) 


and  thus 


S 


L'^Mdc 


(10) 


is  the  scattering  matrix  of  the  problem  in  Figure  1.  While  this  is  quite  true  in 
a  formal  sense,  this  course  leads  to  numerical  difficulties.  The  scattering 
matrix  in  Equation  10  exhibits  the  relative  convergence  phenomenon  (see  the 
Introduction  and  Brief  History  section)  and  is  unsuitable  for  direct  numerical 
solution  of  the  unknown  outputs. 

By  satisfying  the  boundary  conditions  in  a  least-squares  sense,  one  avoids 
the  relative  convergence  problem.  In  terms  of  the  scattering  matrix,  we  use 
the  fact  that  the  adjoint  (complex  conjugate  transpose)  of  the  L  matrix  can  be 
multiplied  through  from  the  left  on  each  side  of  Equation  8.  Thus,  from 
Equation  8, 


L‘  L  b  =  U  M  a  (11) 

Integrating  L^L  and  L^M  over  the  appropriate  parts  of  the  boundary  at  z  =  0 
results  in 


S  =  H-lf 


(12) 


where  H  =  Jc  L^L  dc  is  square  and  Hermitian,  and  f  =  Jc  L‘M  dc  is  the  known 
integrated  input  (sometimes  called  the  forcing  function).  Of  course,  the 

scattering  matrix  of  Equation  12  reduces  formally  to  that  of  Equation  10,  but 
Equation  12  is  numerically  stable.  The  scattering  matrix  of  Equation  12  is 
guaranteed  to  converge  to  a  given  result  as  the  number  of  modes  in  either 
region  is  increased.  The  scattering  matrix  in  Equation  10  will  converge  to 
different  answers  depending  upon  the  number  of  modes  used  and  the  modal 
ratio  r  =  m/p  as  m  and  p  arc  increased  (Reference  4). 
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The  error  associated  with  this  least-squares  scattering  matrix  is  found 
from  Equation  1.  The  "distance"  between  two  complex  quantities  in  a  physical 
sense  (References  16  and  17)  is  determined  by 


d(f,g)  =  lf-gl  (13) 

Thus, 

e  =  lb-Sal  (14) 

gives  a  physically  meaningful  indication  of  how  far  away  (i.e.,  when  m  and  p 
approach  infinity)  any  truncated  result  is  from  the  exact  answer.  Thus, 


e  =  (b^  -  a‘S‘)  (b  -  Sa) 


=  b‘b  -  b^Sa  -  a^S^b  +  a‘S‘Sa 


(15) 


If  the  boundary  conditions  arc  perfectly  matched,  the  error  e  is  zero.  It  may 

not  be  possible  to  have  c  =  0  exactly  except  when  the  number  of  modes  taken  on 

either  side  of  the  junction  is  allowed  to  approach  infinity.  Since  truncation  is 

a  numerical  necessity,  the  set  of  unknown  outputs  [i.e.,  bm^*^  and  bp(2)]  that 
causes  the  error  to  be  a  minimum  is  the  physically  correct  set  (References  18 

and  19). 


CANONICAL  PROBLEMS 


We  depart  from  the  theory  of  the  General  Theory  section  to  discuss  two 

particular  waveguide  discontinuities:  (1)  the  infinitely  thin  capacitive  iris  in 
a  parallel-plate  waveguide  and  (2)  the  asymmetrical  H-plane  parallel-plate 
waveguide  step.  The  above  were  chosen  as  test  cases  for  a  number  of  reasons; 
(a)  Both  arc  two-dimensional — this  means  there  will  be  fewer  nonzero  field 
components  to  match  across  the  discontinuities;  (b)  both  have  been  done 

before  by  other  methods  so  that  comparison  values  exist;  and  (c)  (1)  and  (2) 
are  in  different  problem  classes  and  thus  provide  a  fair  test  of  the  LSBRM. 

The  details  of  deriving  the  H  matrices  that  must  be  inverted  to  compute  the 

reflection  and  transmission  coefficients  for  (1)  and  (2)  arc  given  in 

Appendixes  A  and  B,  respectively. 
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In  the  capacitive  iris  problem  (see  Figure  A-1),  we  have  a  number  of 
conditions  that  ensure  the  validity  of  the  solution.  First  of  all,  each  of  the 
tangential  boundary  conditions  must  be  satisfied  by  the  computed  coefficients, 
i.e..  Equations  A-7  must  be  satisfied  simultaneously.  Thus,  the  tangential 
electric  fields  on  either  side  of  the  iris  from  b  <  y  <  a  must  go  to  0  and  be 
continuous  from  0  <  y  <  b. 

Consider  Figure  2,  where  the  waveguide  height,  a,  is  1.2  centimeters  and 
the  iris  has  a  length  b  =  a/2  =  0.6  centimeter  at  a  frequency  of  10  gigahertz 
(Xo  =  3  centimeters).  This  is  a  plot  of  the  normalized  tangential  electric  field 
magnitudes  in  regions  1  and  2  [i.e.  lEy(l)|and  lEy(2)|]  as  functions  of  the  guide 
height  along  the  discontinuity  plane  at  z  =  0.  It  is  clear  that  using  26  modes  in 
each  region,  i.e.,  M  =  N  =  26,  gives  virtually  identical  fields  across  the  aperture 
and  causes  the  fields  to  drop  sharply  to  0  at  y  =  b  =  0.6  on  the  perfectly 

conducting  iris.  The  real  and  imaginary  parts  of  the  dominant  transverse 

electromagnetic  (TEM)  mode  reflection  and  transmission  coefficients  are 

given  as  well  as  the  normalized  susceptance  B/Yq  of  the  associated  equivalent 
circuit.  The  two  problems  with  Figure  2  are  that  the  fields  are  not  0  at  exactly 
y  =  0.6,  and  B/Yq  =  1.725,  which  is  in  error  by  8.5%.  By  taking  a  larger  number 
of  modes,  we  can  correct  these  problems.  For  76  modes  on  each  side  with  all 
other  parameters  the  same  as  in  Figure  2,  the  fields  fall  to  exactly  0  at  y  =  0.6 
and  the  normalized  susceptance  is  in  error  by  only  3.8%,  as  shown  in  Figure  3. 

Similarly,  the  tangential  magnetic  fields  must  be  continuous  in  the 
aperture  and  differ  from  each  other  by  the  current  density  on  the  iris. 

Looking  at  Figures  4  and  5.  we  see  that  this  is  true.  Figure  4  compares  the 

tangential  magnetic-field  magnitudes  for  26  modes  on  each  side  of  the  iris.  In 

Figure  5  with  76  modes  in  each  region,  the  fields  match  better,  particularly  at 
the  edge  of  the  iris  (y  =  0.6),  and  the  fields  in  the  aperture  are  tending  toward  0 
as  the  number  of  modes  increases.  By  symmetry,  we  know  that 
=  0  for  0  <  y  <  b  if  the  boundary  conditions  are  matched  exactly.  Thus, 

for  increasing  mode  number,  we  note  that  it  is  more  difficult  to  match  the 

tangential  magnetic  fields  over  the  discontinuity  plane  than  the  tangential 

electric  fields,  and  further  experimentation  shows  that  this  is  always  true. 
This  may  be  a  consequence  of  the  fact  that  the  tangential  boundary  conditions 
on  the  magnetic  fields  o’^er  the  perfectly  conducting  portion  of  the 
discontinuity  [i.e.,  h  x  -  H^^))  =  Jj]  arc  not  explicitly  used  in  the  original 

formulation  of  the  scattering  matrix. 

Now  that  the  tangential  boundary  conditions  have  been  satisfied,  we  must 
check  that  normal  conditions  are  also  satisfied,  although  this  is  usually  taken 
for  granted.  In  this  problem,  since  the  iris  generates  only  higher  order 

transverse  magnetic  (TM)  modes,  the  only  nonzero  normal  field  component  to 
the  boundary  at  z  =  0  is  Ez-  We  know  that  continuous  normal  B  and  D  arc  the 

normal  conditions  to  be  satisfied  across  the  aperture;  therefore,  on  0  <  y  <  b, 

Ej(^)  =  Ez^^).  Looking  at  Figures  6  and  7,  we  see  that  this  is  indeed  true  for  26 

and  76  modes  in  each  region,  respectively.  In  fact,  using  V  •  D  =  0  and  the 
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continuity  equation,  for  0  <  y  <  a  and,  thus,  continuous  normal  D  is 

true  along  the  entire  discontinuity  plane. 

The  third  condition  to  be  met  is  that  for  a  =  0.4  A,  and  b  =  a/2,  the  normalized 
susceptance  across  the  junction  is  B/Yo  =  1.59.  From  Figure  8,  we  see  that  as 
the  number  of  modes  is  increased,  the  susceptance  is  tending  toward  this  value 
for  a  modal  ratio,  M/N  =1.  At  75  modes  on  each  side,  the  susceptance  value  is 
1.65,  which  is  -  3.8%  error  as  compared  with  the  exact  value. 

As  we  mentioned  in  the  General  Theory  section,  the  error  function  can  be 
computed  from  Equation  15  when  the  coefficients  are  known.  This  error  must 
be  a  minimum  for  a  given  set  of  modes,  and  we  want  it  to  tend  to  0  as  the 
number  of  modes  increases.  From  Figure  9  for  a  modal  ratio  of  1,  we  see  that 
the  error  function  for  the  reflection  and  transmission  coefficients  (i.e.,  the 
b's)  is  less  than  5%  after  only  8  or  10  modes  on  each  side,  and  it  decreases 
toward  0  in  a  smooth,  well-behaved  manner. 

In  order  to  compare  the  theoretical  model  with  experiment,  the  quantities 
of  interest  are  the  magnitude  and  phase  of  the  dominant  mode  reflection  and 
transmission  coefficients.  Figure  10  shows  the  magnitude  of  the  reflection 
and  transmission  coefficients  as  functions  of  increasing  mode  number  for  a 
modal  ratio  of  1.  Figure  11  shows  the  phases  as  functions  of  increasing  mode 
number.  After  20  or  25  modes,  both  quantities  are  changing  only  slightly  with 
increasing  mode  number.  These  two  plots  confirm  the  results  of  the  error 

function  in  Figure  9.  For  practical  purposes,  it  is  necessary  to  take  only  10  or 
15  modes  on  each  side  of  the  discontinuity  to  obtain  less  than  5%  accuracy  in 
the  dominant  mode  reflection  and  transmission  coefficient  values.  It  is  also  of 
interest  to  compute  the  reflection  and  transmission  coefficients  over  a 
frequency  band.  Figures  12  and  13  show  the  reflection  and  transmission 
coefficient  magnitudes  and  phases,  respectively,  at  X-band.  Over  the 
frequency  range  shown,  the  reflection  coefficient  magnitude  increases  with 
increase  in  frequency,  whereas  the  transmission  magnitude  decreases.  The 
phases  of  both  coefficients  remain  relatively  constant  over  this  range. 

Finally,  although  we  have  shown  that  the  numerical  results  all  converge 
as  the  number  of  modes  increases,  we  must  also  show  convergence  with 
respect  to  the  modal  ratio.  So  far  the  modal  ratios  in  all  cases  have  been  1.  But 
consider  how  the  normalized  susceptance  is  affected  by  the  modal  ratio.  In 
Figure  8,  where  the  ratio  is  1  at  25  modes  on  each  side,  the  susceptance  is  about 
1.725.  In  Figure  14,  where  the  ratio  is  2  for  50  modes  in  region  1  and  25  modes 
in  region  2,  B/Yq  =  1.735.  In  Figure  15,  where  the  ratio  is  3  for  75  modes  in 
region  1  and  25  modes  »"  region  2,  B/Yq  =  1.74.  Obviously,  while  the 

susceptance  is  tending  to  the  correct  answer  regardless  of  modal  ratio,  a 
poorly  chosen  ratio  will  cause  a  much  larger  number  of  total  modes  to  be  used 
in  order  to  achieve  a  particular  accuracy.  In  this  problem,  since  the 

waveguide  is  the  same  height  on  either  side  on  the  iris,  intuition  indicates  that 
an  equal  number  of  modes  in  each  region  is  the  "best"  ratio.  More  will  be  said 
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concerning  this  aspect  of  the  analysis  for  the  H-plane  parallel  plate  step 
problem. 

Thus,  after  considering  all  of  the  above  checks  on  the  capacitive  iris 
solution,  we  can  conclude  that  the  LSBRM  is  numerically  stable,  accurate  to 
less  than  5%  (depending  on  the  number  of  modes  one  wishes  to  take),  and  is 
slowly  convergent,  although  the  convergence  appears  to  be  guaranteed. 

In  the  asymmetrical  H-plane  parallel-plate  waveguide  problem  (details  in 
Appendix  B),  the  situation  is  a  little  more  complex  (see  Figure  B-1).  Again,  we 


(at  z  =  0)  have  our  aperture  between  0  <  x  <  b  and  a  perfectly  conducting 
surface  between  b  <  x  <  a.  For  TEqi  mode  incidence,  the  step  generates  only 
higher  order  TE  modes  and,  thus,  Hz,  Hx,  and  Ey  are  the  only  nonzero  field 
components.  The  tangential  boundary  conditions  on  z  =  0  that  arc  used  in 
formulating  the  H-matrix  are 

=  0  ;  b  <  x  <  a 

(16a) 

;  0  <  X  <  b 

(16b) 

=  ;  0<x<b 

(16c) 

An  auxiliary  condition  not  explicitly  included  is  the 
condition  on  the  step.  (For  further  discussion  of 

Appendix  C.) 

tangential  magnetic  field 
boundary  conditions,  see 

u(l)  Tjd)  T  .  U  ^  ^ 

Hx  =H,an  =  Js  >  b<x<a 

(17) 

[We  assume  that  and  £(2)  no  longer  exist  in  region  2  when  x  >  b,  since 

that  region  is  bounded  by  perfect  conductors.] 

The  normal  conditions  on  z  =  0  are 


(1)  _  ti(2) 

z  “z 

;  0  <  X  <  b 

(18a) 

(1) 

z  =  0  ; 

b  <  X  <  a 

(18b) 

Not  only  must  all  of  the  above  be  satisfied,  but  the  error  function  must  be  a 
minimum  as  before. 
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To  show  that  the  previous  conditions  are  all  satisfied,  we  use  Figure  B-1 
and  let  (1)  the  larger  guide  height,  a,  be  2.7  centimeters,  (2)  the  smaller  guide 
height,  b,  be  2  centimeters,  and  (3)  the  frequency  be  10  gigahertz.  M  and  N 
refer  to  the  maximum  number  of  modes  used  in  the  large  and  small  guides, 
respectively.  Figures  16  and  17  are  plots  of  the  electric  field  magnitudes  as 
functions  of  waveguide  height  along  the  step  at  z  =  0.  Figure  16  shows  how 
well  the  fields  match  for  20  modes  in  each  region,  while  Figure  17  shows  the 
match  for  50  modes  in  each  region.  For  50  modes,  the  tangential  electric  field 
boundary  condition  is  satisfied  almost  exactly. 

The  tangential  magnetic  field  magnitudes  as  functions  of  the  guide  height 
along  z  =  0  are  shown  in  Figures  18  and  19.  As  for  the  capacitive  iris  problem, 
the  tangential  magnetic  fields  are  more  difficult  to  match  than  the  tangential 
electric  fields.  In  Figure  18,  lHx^'^1  and  are  fairly  well  matched  for  20 

modes  until  close  to  the  discontinuity,  where  IH  x^^^  I  goes  to  0  (it  no  longer 

exists),  but  Hx^*^=  Js-  Figure  19,  the  fields  match  better  for  50  modes  in 
each  region;  however,  there  is  still  some  discrepancy  close  to  the 
discontinuity. 

The  only  normal  condition  is  continuous  normal  B  (or  H,  in  this  case)  on 

0  <  X  <  b,  and  =  0,  as  shown  above,  on  b  <  x  <  a.  Figures  20  and  21  are  plots 

of  the  normal  magnetic  field  magnitudes  as  functions  of  the  guide  height  at  z  = 
0.  For  20  modes  in  each  region  (Figure  20),  the  field  match  is  poor  in  both  the 
aperture  region  and  the  step  region.  Not  until  we  increase  the  number  of 
modes  in  each  region  to  75  (Figure  21),  do  the  normal  fields  match  well. 

Figure  22  is  a  plot  of  the  error  function  versus  mode  number  for  a  modal 
ratio  of  1,  Since  the  step  in  this  example  is  not  very  large,  it  is  necessary  to 
take  only  7  or  8  modes  to  obtain  an  accuracy  of  less  than  4%  in  the  computed 
reflection  and  transmission  coefficients.  The  dominant  mode  reflection 
coefficient  magnitude  as  a  function  of  mode  number  is  shown  in  Figure  23. 
This  plot  supports  the  results  of  the  error  function,  i.e.,  only  7  or  8  modes  need 
to  be  included  to  give  accurate  values  of  reflection. 

Now  consider  the  example  where  a  =  2.7  but  b~  a/2  =  1.4  at  10  gigahertz. 
This  is  a  much  larger  step  than  in  the  previous  example,  and  it  will  serve  to 
illustrate  an  interesting  effect  concerning  the  modal  ratio.  Consider  Figures 
24  through  26.  If  the  total  number  of  modes,  P  =  M  +  N,  is  100,  then  Figures  24, 
25,  and  26  show  the  error  as  a  function  of  increasing  mode  number  for  modal 
ratios  1,  2,  and  3,  respectively  (i.e.,  M  =  50,  N  =  50;  M  =  66,  N  =  33;  and  M  =  75,  N  = 
25).  The  geometric  ratio  for  the  waveguide  heights  is  a/b  ~  2.  We  would  like  to 
find  the  smallest  possible  error  for  a  given  value  of  P.  By  comparing  Figures 
24  through  26,  it  is  clear  that  a  modal  ratio  of  2  gives  the  smallest  error  when 
P  =  100.  Also,  we  would  like  to  find  the  smallest  possible  error  for  the  smallest 
number  of  total  modes  P.  For  example,  M  =  50,  N  =  50  gives  an  error  of  3.83%, 
whereas  M  =  50,  N  =  25  gives  an  error  of  3.96%.  An  increase  in  accuracy  of 
0.13%  is  not  worth  the  extra  computation  needed  to  include  25  extra  modes.  We 
see  that  the  modal  ratio  is  controlled  in  some  sense  by  the  geometric  ratio,  a/b. 


13 


NWC  TP  6941 


and  in  general  we  can  confirm  the  known  result  that  if  a/b  ~  2,  then  rn/n  ~  2  is 
the  "best"  modal  ratio  both  in  terms  of  the  error  and  the  total  number  of  modes 
P. 


CONCLUSIONS 


We  have  used  the  least-squares  boundary  residual  method  (LSBRM)  to 
obtain  very  accurate  numerical  values  for  the  modal  reflection  and 
transmission  coefficients,  the  electric  and  magnetic  field  components,  and  the 
error  associated  with  certain  types  of  waveguide  discontinuity  problems.  Wc 
have  considered  the  general  theory  in  terms  of  a  scattering  matrix 
representation  and  have  shown  that  the  theory  works  very  well  for  two 
canonical  problems:  (1)  the  thin  capacitive  iris  in  a  parallel-plate  waveguide 
and  (2)  the  asymmetrical  H-plane  parallel-plate  waveguide  step.  The  method 
exhibits  slow  convergence  for  large  discontinutics;  however,  its  advantages — 
explicit  evanescent  modes,  guaranteed  convergence,  excellent  accuracy,  and 
an  associated  error  function,  which  must  be  a  minimum  for  the  physically 
correct  solution — arc  significant. 

In  the  immediate  future,  we  will  continue  to  look  at  similar  discontinuity 
problems  but  with  the  following  modifications: 

1.  Discontinuities  in  cavities  with  perfectly  conducting  end  plates 

2.  Discontinuities  in  cylindrical  and  spherical  geometries 

3.  Discontinuities  that  generate  both  transverse  electrie  and  transverse 
magnetic  higher  order  modes 

4.  Multiple  discontinuities  less  than  a  wavelength  apart 

This  is  a  very  powerful  and  versatile  technique  for  predicting  the 
scattering  properties  of  any  discontinuity  or  junction,  provided  the  geometry 
can  be  represented  in  terms  of  cigcnscrics.  It  will  enhance  our  capability  to 
model  complex  microwave  devices  so  that  their  scattering  properties  can  be 
well  understood. 
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FIGURE  1.  Geometry  for  Junction  Between  Two  General  Waveguides. 
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FIGURE  2.  Normalized  Tangential  Electric  Field  Magnitudes  Versus  Guide 
Height  for  Capacitive  Iris  Junction  (26  Modes). 
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FIGURE  3.  Normalized  Tangential  Electric  Field  Magnitudes  Versus  Guide 
Height  for  Capacitive  Iris  Junction  (76  Modes). 
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FIGURE  4.  Normalized  Tangential  Magnetic  Field  Magnitudes  Versus  Guide 
Height  for  Capacitive  Iris  Junction  (26  Modes). 
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FIGURE  5.  Normalized  Tangential  Magnetic  Field  Magnitudes  Versus  Guide 
Height  for  Capacitive  Iris  Junction  (76  Modes). 
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FIGURE  6.  Normalized  Normal  Electric  Field  Magnitudes  Versus  Guide  Height 
for  Capacitive  Iris  Junction  (26  Modes). 
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FIGURE  7.  Normalized  Normal  Electric  Field  Magnitudes  Versus  Guide  Height 
for  Capacitive  Iris  Junction  (76  Modes). 
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FIGURE  8.  Normalized  Susccptance  Versus  Mode  Number  for  Capacitive  Iris 
Junction  (Modal  Ratio  of  1). 
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FIGURE  9.  Error  Versus  Mode  Number  for  Capacitive  Iris  Junction  (Modal  Ratio 
of  1). 
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FIGURE  10.  TEM  Mode  Reflection  and  Transmission  Magnitude  Versus  Mode 
Number  for  Capacitive  Iris  Junction. 
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FIGURE  11.  TEM  Mode  Reflection  and  Transmission  Phase  Versus  Mode  Number 
for  Capacitive  Iris  Junction. 
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FIGURE  12.  TEM  Mode  Reflection  and  Transmission  Magnitude  Versus 
Frequency  for  Capacitive  Iris  Junction. 


20 


NWC  TP  6941 


8.0  8.5  9.0  9.5  10.0  10.5  11.0  11.5  12.0 


FREQUENCY,  GHz 

FIGURE  13.  TEM  Mode  Reflection  and  Transmission  Phase  Versus  Frequency 
for  Capacitive  Iris  Junction. 
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FIGURE  14.  Normalized  Susccptancc  Versus  Mode  Number  for  Capacitive  Iris 
Junction  (Modal  Ratio  of  2). 
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FIGURE  15.  Normalized  Susceptance  Versus  Mode  Number  for  Capacitive  Iris 
Junction  (Modal  Ratio  of  3). 
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FIGURE  16.  Normalized  Tangential  Electric  Field  Magnitudes  Versus  Guide 
Height  for  Parallel-Plate  Waveguide  Junction  (TEqi  Mode  Incidence  From 
Larger  Guide,  20  Modes). 
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FIGURE  17.  Normalized  Tangential  Electric  Field  Magnitudes  Versus  Guide 
Height  for  Parallel-Plate  Waveguide  Junction  (TEqi  Mode  Incidence  From 
Larger  Guide,  50  Modes). 
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FIGURE  18.  Normalized  Tangential  Magnetic  Field  Magnitudes  Versus  Guide 
Height  for  Parallel-Plate  Waveguide  Junction  (TEqi  Mode  Incidence  From 
Larger  Guide,  20  Modes). 
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FIGURE  19.  Normalized  Tangential  Magnetic  Field  Magnitudes  Versus  Guide 
Height  for  Parallel-Plate  Waveguide  Junction  (TEqi  Mode  Incidence  From 
Larger  Guide,  50  Modes). 


GUIDE  HEIGHT,  cm 
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= 

h;2), 
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= 

20 

A 
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= 
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= 

0.883 

T. 

1 

= 

-0.014 

FIGURE  20.  Normalized  Normal  Magnetic  Field  Magnitudes  Versus  Guide 
Height  for  Parallel-Plate  Waveguide  Junction  (TEqi  Mode  Incidence  From 
Larger  Guide,  20  Modes). 
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GUIDE  HEIGHT,  cm 


a 

• 

= 

h'(2) 

z 

M 

75 

N 

= 

75 

A 

= 

2.70  cm 

B 

= 

2.00  cm 

FREQ 

= 

10.0  GHz 

R 

-0.070 

r 

Ri 

= 

-0.121 

T 

_ 

0.883 

r 

T. 

= 

-0.001 

FIGURE  21.  Normalized  Normal  Magnetic  Field  Magnitudes  Versus  Guide  Height 
for  Parallel-Plate  Waveguide  Junction  (TEqi  Mode  Incidence  From  Larger 
Guide,  75  Modes). 
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B 

FREQ 
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2.70  cm 
2.00  cm 
10.0  GHz 


FIGURE  22.  Error  Versus  Mode  Number  for  Parallel-Plate  Waveguide  Junction 
(TEqi  Mode  Incidence,  Modal  Ratio  of  1). 
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MODE  NUMBER,  N 


FIGURE  23.  Reflection  TEoi  Mode  Coefficient  Magnitude  Versus  Mode  Number 
for  Parallel-Plate  Waveguide  Junction. 


MODE  NUMBER.  N 


FIGURE  24.  Error  Versus  Mode  Number  for  Parallel-Plate  Waveguide  Junction 
(Modal  Ratio  of  1). 
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MODE  NUMBER.  N 


FIGURE  25.  Error  Versus  Mode  Number  for  Parallel-Plate  Waveguide  Junction 
(Modal  Ratio  of  2). 


MODE  NUMBER.  N 


FIGURE  26.  Error  Versus  Mode  Number  for  Parallel-Plate  Waveguide  Junction 
(Modal  Ratio  of  3). 


NWC  TP  6941 


Appendix  A 

A  CAPACITIVE  IRIS  IN  A  PARALLEL-PLATE  WAVEGUIDE 


Consider  the  simple  capacitive  iris  discontinuity  in  a  symmetric  parallel- 
plate  waveguide  (References  13  and  20)  (see  Figure  A-1). 


— b-TEM 
REGION  1 


REGION  2 


FIGURE  A-1.  Geometry  for  the  Infinitely  Thin  Capacitive  Iris  in  a 
Parallel-Plate  Waveguide. 


The  following  assumptions  are  made:  (1)  Only  transverse  electromagnetic 
(TEM)  mode  incidence  exists  (from  the  left);  (2)  both  regions  extend  to 
infinity;  (3)  only  higher  order  transverse  magnetic  (TM)  modes  arc  generated 
by  the  discontinuity;  (4)  all  3/3  x  =  0,  since  the  waveguide  extends  to  infinity 
along  the  x  coordinate;  and  (5)  e  =  Eq  and  =  |Uo  inside  the  guide. 

Assumption  (1)  indicates  that 


a;^^  =  0  ;  m>0 


(A-1) 


(The  value  of  1  is  arbitrarily  chosen.) 
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Assumption  (2)  implies  that 


Vp 


(A-2) 


The  only  nonzero  field  comp  tnents  are  E^,  Hx,  and  Ey  from  assumptions  (3)  and 
(4).  The  normalized  tangential  fields  in  region  1  are  well  known  to  be 
(Reference  13) 


a  ao  e  + 


m=0 


mtcy 

cos - e'™ 

a 


(A-3a) 


(A-3b) 


while  the  fields  in  region  2  are 


pity 

cos - e  'p 

a 


(A-4a) 


ZoH 


(2) 

X 


p=0 


Yf  cos 


p^ye^^ 

a 


z>  0 


(A-4b) 
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where 


1  ;  j  =  0 


if k  =  l,j  =  morif  k  =  2, j  =  p 


2  ;  j  >  0 


(A-5) 


and  the  normalized  wave  admittances  in  each  region  are 


;  ifk=  l,j  =  morifk  =  2,j  =  p 


(A-6) 


(Zq  is  the  frce-space  impedance.)  The  boundary  conditions  at  z  =  0  arc 


=  0  ;  b<y<a 

(A-7a) 

^^^  =  0  ;  b<y<a 

(A-7b) 

=  ;  0<y<b 

(A-7c) 

=  H®  ;  0  <  y  <  b 

(A-7d) 
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Substituting  Equations  3  and  4  into  7,  the  boundary  conditions  become 


Let 


y  <  b 


m=0 


o(2) 

^?\/^cos-^  =  0  ;  b<y<a 


p=0 


^/4  =<  ■  \/4-=n® 

cil>  =  cosS  ,  = 


0  <  y  <  b 


(A-8a) 


(A-8b) 


(A-8c) 


(A-8d) 


(A-9a) 


(A-9b) 
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In  matrix  form.  Equation  8  becomes 


'  N«ci:> 

T^(l)  ^(l) 


N, 


(1)  p(i) 
m  ''Tn 

0 


_n(2)c(2)  ^ 


N 


(2)  p(2)  y(2) 
0 


\  rr* 

r-n 

bm 

/i  ,(i) 

1 

VDp  J 

"V 

-1 

1  0. 

o  r 


Lb  =  Ma 


The  adjoint  matrix  to  L  is 


L+  = 


1^(1)  p(l)  vt(1)  p(l)  0 


^(2)  p(2)  y^2)^ 

i>p'  '-p-  Ip' 


0 


N 


(2)  ^(2) 


■  / 


(A-IO) 


(A-11) 


(m',  p'  are  dummy  indices). 

The  product  L^L  must  be  integrated  over  the  appropriate  portion  of  the 
cross-sectional  boundary  where  each  of  the  original  boundary  conditions 
holds.  Let 


H  =  J^L^Ldc 


(A-12) 


where  C  is  the  entire  boundary  at  z  =  0.  We  can  write  H  in  partitioned  form  as 


H=i 


VH 


(21) 


H 


(12)^ 


h(22)^ 


(A-13) 
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is  formed  by  multiplying  the  first  row  of  by  the  first  column  of  L  and 
integrating,  noting  that  m  =  0  to  «  ;  m'  =  0  to  oo. 


h'”> 


Jo  Jo 

+  r  n“.’  n“>  ci;.>  c2’  dy 
Jb 


^(1)  p(l)  p(l) 

^m’ 


dy 


(A-14) 


h(ii) 


(A-15) 


(See  the  end  of  this  Appendix  for  evaluation  of  integrals.)  Thus, 


_  If  1  T 

“  ~  *Sn  *ni  *m'  •'mm' 


Without  detail. 


H“‘’  =  V«.tY“’'Yi;’-Il 


H“"=Kp+ y'?’‘ Y®Jpp. 


(A-16a) 


(A-16b) 

(A-16c) 

(A-16d) 


Thus  and  H(22)  are  matrices  such  that  if  m.  m'  run  from  0  to 

M  and  p,  p'  run  from  0  to  P,  then  has  dimensions  (M  +  1)  by  (M  +  1), 

is  (M  +  1)  by  (P  +  1),  H(21)  is  (P  +  1)  by  (M  +  1),  and  H(22)  is  (p  +  i)  by  (P  +  1). 

Of  course,  when  M  =  P, 


H^l)  =  H<22) 


(A-I7a) 
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and 


H(12)  =  h(21)  (A-17b) 

because  of  the  particular  symmetry  of  this  problem.  To  integrate  the  forcing 
function,  we  take 


•=Jl* 


Mdc 


which  is  a  column  vector  (sec  end  of  this  Appendix  for  integrals), 


f  = 


with 

p^. 


(A-18) 


(A-19a) 

(A-19b) 


Once  H  and  f  arc  known,  we  can  take  the  inverse  of  H  and  find  the  scattering 
matrix.  Since  only  a  single  input  was  assumed,  we  obtain  only  one  row  and 
one  column  of  the  complete  scattering  matrix; 


S  =  H'l  f 


(A-20) 


From  previously,  if 


H  b  =  f  a 


(A-21) 
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the  set  of  unknown  coefficients,  and  can  be  found  independently 

using  the  partitioned  form  of  H  (References  21  and  22).  Thus, 


(  jjdi)  ’  Y  Ji)  ^ 


\P  ) 


(  j(i)  ^ 
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(Let  =  1-)  Then 


(A-23a) 


^(21)  jj(l)  ^  j^(22)  jj(2)  _  j(2) 


(A-23b) 


and 


5(2)  _  Jj^(22)  _  J^(21)  H^12)j-1  JJ<2)_  j^(21) 


(A-24a) 

(A-24b) 


The  reflection  and  transmission  coefficients  are  b(')  and  b^^),  respectively,  for 
the  TEM  plus  all  higher  order  TM  modes.  For  computing  purposes,  rather  than 
use  Equations  24,  it  is  easier  to  use 


b  =  H-l  f 


(A-25) 


and  find  b^^^  and  b^^^  simultaneously.  H  is  an  (M  +  P  +  2)  by  (M  +  P  +  2)  matrix, 
while  f  is  an  (M  +  P  +  2)  column  vector  as  is  b. 
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ORTHONORMALIZED  INTEGRALS  USED  IN  THE  CAPACITIVE  IRIS  PROBLEM 


Let 


■^mm'  ~ 


-f 


m  Tcy  mTty  . 

cos - —  cos - dy 

a  & 


(A-26) 


JiniTi’ 


b 

a 


;  m  =  0,  m'  =  0 


VT  .  mTtb  ,  ^ 

- sin -  ;  m  =  0,  m  0 

mjc  a 

VT  .  m'Ttb  .  -  _ 

— sin - ;  m  ^  0,  m  =  0 

m  7C  a 


;  m  =  m'  (m,  m'  >  0) 


sin 

f  m'K 

mre^ 

lb  .infiSil+iSi'lb 

1. 

1  a 

a  ^ 

|b  +  — jb 

a 

“7 

m'TC 

mn) 

m'Tc  mn 

l  V 

.  a 

a  J 

a  a 

m  ^m' ; 


(A-27) 


V  f  m'Tty  mny 


I 


cos 


cos - ^  dy  = 

a  a 


fO  m?tm' 

1  m  =  m', 
m,  m'  >  0 
1  m  =  m'  =  0 


(A-28) 


•^m'p 


f 

Jo 


m  Ky  prey  , 

cos - -  cos  — -  dy 

a  a 


(A-29) 
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(Same  as  Equation  A-27  for  Jmm'  with  m'  =  m,  p  =  m'.) 

f 

Jo 


^p’m  ~ 


mTCy  p  jcy  , 
cos  — -  cos-^^ — —  dy 
0  a  a 


(A-30) 


(Same  as  Jmm'  with  p'  =  m.  m  =  m'.)  Jpp*  =  Jmm'  with  p  =  m,  p'  =  m' 


o  r 


0  ;  m'  >  0 


1  ;  m’  =  0 


Pm'  = 


cb 

^m'  f  mjty  , 

cos  — -  dy 
a  Jo  " 


f  b 
a 


Pm' 


;  m'  =  0 


VT  .  m'Kb  ,  ^ 
-T-  sin  — - —  ;  m  >  0 
(  in  K  a 


(pp’  is  the  same  with  m'  =  p'.) 


(A-31) 


(A-32) 


(A-33) 
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Appendix  B 

THE  ASYMMETRICAL  H-PLANE  PARALLEL-PLATE  WAVEGUIDE  STEP 


Assume  that  the  asymmetrical  H-planc  parallel-plate  waveguide  step  is  as 
shown  in  Figure  B-l.  The  following  assumptions  arc  made;  (1)  There  is  TEq  i 


z  =  0 


FIGURE  B-l.  Geometry  for  the  Asymmetrical  H-Planc  Parallel-Plate  Waveguide 
Step. 


mode  incidence  from  the  left;  (2)  both  regions  extend  to  infinity;  (3)  only 
higher  order  TE  modes  arc  generated  by  the  step;  (4)  all  3/5y  =  0;  and  (5)  e  =  Eq. 
p  =  Po  inside  the  guides. 

We  have  (similar  to  Appendix  A), 


a<”.l 

=  0  ;  m  >  1 


and 


;  Vp 


(B-la) 


(B-lb) 
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Nonzero  components  of  the  fields  are  Hz,  Hx,  and  Ey .  The 
tangential  fields  in  region  1  are 


e  'm 


]  sin 


mKx 

a 


mTtx 

a 


z  <  0 


and  in  region  2,  they  are 


e  p 


(2) 

Y  z 


ptex 

sin 

0 


z  >  0 


7  14(2)  _  n  \  Y(2)h<2),.-Yp  z  P’TX 

P=1 

with  the  axial  propagation  constants  given  by 


and  the  normalized  admittances  are 


normalized 


(B-2a) 


(B-2b) 


(B-3a) 


(B-3b) 


(B-4) 
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The  tangential  boundary  conditions  used  on  z  =  0  arc 
4^^  =  0  ;  b  <  X  <  a 

=  ;  0<x<b 

=  ;  0<x<b 

Substituting  Equations  B-2  and  B-3  into  B-6,  we  have 


p=i 


(B-5) 

(B-6a) 

(B-6b) 

(B-6c) 

(B-7a) 
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sin 


Ttx  ^(1)  .  mine 

-r  +  2^  b;„  sm 


m=l 


=  0 


(B-7c) 


In  matrix  form. 


L  = 


Vi“ 

VI 


ItlTCX 


Fl  .  pTlX 

“V  b  b 

V  b  ^  b 


/T  .  mjtx 

^  ~  sin - 

V  a  a 


0 


and 


L'  = 


fl  .  m  7CX  /T  ,,(!)*  .  m  Jtx  /T 

*  /  —  sin -  A  /  “  Y' .  sm -  ^  “ 

Va  a  va*"  ava 


sin 


m  Jix 


fl  .  p'TtX  fl  (2)*  .  p' 

V  b  b  V  b  ^P’  i 


7tX 

b 


0 


where  m'  and  p'  are  dummy  indices.  Lei  the  matrix 


(B-8) 


(B-9) 


b  = 


(B-10) 


contain  the  reflection  and  transmission  coefficients  for  each  mode  and  let  the 
inputs  be 
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Ma  =  al 


(1) 


[2  .  m 
-*  /  —  sin  — 
V  a  a 

>/! 


^  Y<»  sin 
a  ‘  a 


/T  .  :cx 
-*  /  —  sin  — 
V  a  a 


;  4»=i 


Let 


H=  f  L^Ldc 

Jc 


As  previously,  we  see  that 


=  I  f —  dx + -  y:^v'  y::/ 


fn'TCx  niTCx  ^  2  ^(1).  ^(1) 


^  m  *  m 


4 

Jo 


in  Jtx  .  mnx  , 

sin - sin - dx 

a  a 


1  1  y(1)  ■y(l)*  T  •  m  m' 

1  4-  Yj^  ■'mm  »  m  —  m 


■yCl)*  T 

1.  *m'  *  m  •’mm' 


;  m  m’ 


(For  integrals,  sec  end  of  this  Appendix.)  In  similar  fashion. 


(B-Il) 

(B-12) 

(B-13) 

(B-14a) 

(B-14b) 

(B-14c) 
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h(22)^ 


1  +  Y„  J„„  ;  p  -  p 


‘p  'pp  ’ 


y(2)*y<2), 

Ip'  Ip  Jpp' 


;  p’^p' 


Again,  let 


f=  f  L^Mdc 
Jc 


where 


fyd)  T 
il  Im'  •'Im' 


;  1 


"1  +  Yi  Yj  ;  m  =  l 


=  Tip.  +  1] 


and 


f  = 


vn 


(B-14d) 


(B-15a) 

(B-15b) 

(B-16) 


Once  H  and  f  are  found,  we  can  compute  the  coefficients  for  the  parallel-plate 
step  just  as  in  Appendix  A. 
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ORTHONORMALIZED  INTEGRALS  IN  THE  PARALLEL-PLATE  STEP  PROBLEM 


Let 


„  2  f  .  mjcx  .  micx  , 

Km  =  —  I  sin - sin - dx 

a  Jo  a  a 


1  ;  m  =  m' 

0  ;  m  ?!:  m' 


(B-17) 


,  2  .  m'jTx  .  nwtx  . 

W  =  sin  — sm  — dx 


sin  {m'  -  m)  —  sin  (m'  +  m)  “ 
a  a 


(m'  -  m)  — 
a 


(m'  +  m)  — 
a 


m  m’ 


b 

a 


sin 


2m7tb 


2m7i 


ni  =  m' 


(B-18) 
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2  f”  .  m'TCX  .  pTCx  . 
=  -==r  sin - sin  ox 

•"P  Jo  a  b 


1 

sin 

/  m'jc  pit  > 

la  b  J 

j  b  sin 

f  m'Jt  p7i ' 
1  a  b  > 

m' 

'  m'TC  pjc  ^  / 

.a  b  J  ' 

'  m'Ti  pit  \ 

.  a  b  J 

’  T 

a  sin 


2m'jib 


b- 


2m'jc 


m 

a 


P 

b 


(Tmp'  and  Jpp’  arc  same  wilh  proper  correlation  of  indices.) 


(B-19) 
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Appendix  C 

BOUNDARY  CONDITIONS  FOR  THE  PARALLEL-PLATE  WAVEGUIDE 


In  general,  if  only  a  TEqi  mode  is  allowed  to  propagate  in  Figure  C-1,  there 
are  only  three  nonzero  field  components,  i.e.,  Ey,  Hx,  and  Hz. 


FIGURE  C-I.  Geometry  for  the  Simple  Parallel-Plate  Waveguide  With  No 
Discontinuities. 


If  we  assume  that 


Ey  =  A 


(C-1) 


(since  Ey  must  be  zero  on  both  walls),  then  from  Maxwell’s  curl  equation. 


V 


(C-2) 


H 


X 


1  3Ey 
ico^o  dz 


(C-3a) 


and 


47 


NWC  TP  6941 


1  9E., 

pj  _  ^  Jl 

^  iO)|io  <;>x 

Similarly,  from  the  other  curl  equation, 

p _ LfHi-iSi' 

y  icoeo  L  . 

The  divergence  equations  give 


(C-3b) 


(C-4) 


aEy 

ay 


(C-5a) 


(simply  the  result  of  the  two-dimensional  nature  of  the  problem)  and 


an,  an. 


(C-5b) 


Substituting  Equation  C-1  into  C-3, 


Now  consider  what  happens  when  a  discontinuity  is  introduced  at  z  =  0; 
use 


(C-6a) 


(C-6b) 
i.e.,  we 


48 


NWC  TP  6941 


(a)  continuous  tangential  E  and  H 

(b)  continuous  normal  B  and  D 


for  dielectric  boundaries 


and 

A  A  A  ^ 

(c)  n  X  E  =  0,  n  •  B  =  0,  n  X  H  =  for  perfectly  conducting  boundaries, 

where  Jj  is  the  surface  current  density.  Thus,  at  the  discontinuity  (see 
Figure  B-1), 


(1)  Ey(l)  =  Ey(2) 

(2)  Ey<l)  =  0 

(3)  HxO)  = 

(4)  Hzd)  =  Hz(2) 

(5)  Hjd)  =  0 

(6)  Hx<l)  =  Js 


onz  =  0,  0<x<b 
onz  =  0,  b<x<a 

onz  =  0,  0<x<b 
onz  =  0,  0<x<b 
onz  =  0,  b<x<a 
onz  =  0,  b<x<a 


(C-7) 


(2)  (2) 

(We  assume  that  E  ,  H  do  not  exist  for  x  >  b.)  Conditions  (1),  (2).  (3),  and  (6) 
are  tangential  conditions;  conditions  (4)  and  (5)  arc  normal  conditions  at  the 
discontinuity.  It  is  important  to  realize  that  in  the  least-squares  boundary 
residual  method  (LSBRM),  the  tangential  conditions  arc  used  to  compute  the 
scattering  parameters  and  the  statement  is  often  made  that  the  normal 
conditions  are  "equivalent"  (Reference  14)  to  the  tangential  ones  (the 
implication  being  that  the  normal  conditions  arc  satisfied  automatically). 
While  this  may  be  true  in  some  cases,  in  the  present  instance  we  see  that  there 
is  not  enough  information  from  (4)  and  (5)  alone  to  solve  the  problem. 
Attempts  to  mix  appropriate  tangential  and  normal  conditions  in  the  LSBRM 
formulation  have  not  been  successful. 

In  some  treatments  of  the  parallel-plate  waveguide  step,  the  boundary 
condition 


an. 

- =  0  for  z  =  0,  b  <  X  <  a 

dz 


(C-8) 
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is  used  along  with  the  usual  tangential  conditions  (References  23  and  24).  This 
condition  cannot  be  derived  directly  from  (a),  (b),  or  (c);  however,  indirectly, 
we  recall  from  Equation  C-4  that  if  Ey  =  0  on  some  boundary,  then 


dz  ax 


(C-9) 


also,  implying  that  both  terms  must  either  be  equal  to  each  other  or  both  must 
be  zero  simultaneously.  To  use  3Hx/3z=  0  without  simultaneously  using 
3H2/3x  =  0  is  misleading. 
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